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TECHNICAL NOTE NO. 1692 

DETERMINATION OF BENDING MOMENTS IN PRESSURE -LOADED RINGS 
OF ARBITRARY SHAPE WHEN DEFLECTIONS ARE CONSIDERED 
By F . R . Steinbacher and Hsu Lo 


SUMMARY 


An analytical method has been derived for determining bending- 
moment distribution in rings of arbitrary shape under internal pressure 
loads, with the change of geometric shapes caused by the load being 
considered. For the purpose of clarity, the method developed was 
applied only to double -symmetrical shapes. A differential-integral 
equation has been derived for this purpose and its solution obtained 
in the form of a trigonometric series. 

Charts have been provided for two specific families of rings of 
various proportions and flexibilities. Tests conducted on rings of 
both families agree very well with the analytical calculations. For 
rings belonging to or close to these families, the bending -moment 
distributions and the deflections of the ring can be read directly 
from the curves. For rings of entirely different shapes, an average 
of 20 hours is necessary for the complete solution of the problem. 
Examples have been given to show the method of obtaining the curves. 

On comparing the present results with the results of solutions 
in which deflections have been neglected, it is seen that the bending 
moments previously computed have always been much too conservative. 

The error introduced is considerable when the rings become more and 
more flexible. 

It is believed that by the use of this method curves can be 
drawn for a typical frame in a fuselage and from these curves the 
bending moments and deflections can be estimated for all similar 
frames . 


INTRODUCTION 


Numerous papers have been written regarding the design of the 
rings for monocoque fuselages. All of them were based on the assumption 
that the deflections of the ring' caused by the loads are so small that 
the internal bending moments and shear and axial loads will be unchanged 
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by any small change in the geometric shape of the ring. However, the 
error introduced by neglecting the change in the geometric shape 
becomes more and more important as the size of the airplane increases 
and more flexible rings are used to save weight. It is the object of 
this paper, therefore, to take into account deflections when calcu- 
lating the bending moments in rings. 

The fact that the change of geometric shape has an important 
effect on the final bending-moment distribution in the rings can be 
illustrated by the beam- column analogy of an initially curved bar under 
axial tension loads. For the bar as illustrated in figure 1, the 
initial bending moment is M = Py. If the bar is stiff and deflections 
are small enough to be neglected, the final bending moment is the same . 
But if the bar is flexible and the deflections are comparably large, 
the final bending moment becomes 


M = P(y - 8) 


The difference between these two equations depends on the flexibility 
of the bar and can be considerable for flexible bars. 

In this paper two egressions for the bending-moment distribution 
are derived. The first egression is obtained by considering the 
action of the applied pressure and the final shape of the deflected 
ring. The second expression is obtained from the change of curvature 
of the ring caused by the loads. If these two expressions are set 
equal, a differential- integral equation is obtained. The solution of 
this equation gives the radial deflections of the ring and the bending- 
moment distribution. Another formula is derived from which the angular 
displacements can be determined when the radial deflections are known. 

The rings discussed in this paper are assumed to be under only 
internal pressure. The method recommended herein can be extended to 
any system of external loads and is also applicable to rings of any 
shape although only shapes of double symmetry are discussed in the 
present paper. The following assumptions are made: 

(1) The rings have regular and smooth shapes of the type 
encountered in fuselage frames . 

(2) The thickness of the ring is small compared with its radius. 
Consequently the following formula can be used with an extremely small 
error. (See reference 1.) 


_1_ 

Pl 
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original curvature at certain point of ring 

final curvature at corresponding point of ring after ring 
is deflected 

internal bending moment built up because of deflections 
of ring 

bending flexibility 

( 3 ) The ring is regarded as inextensible . 

(4) The deflections of the ring are defined by its radial 
deflections v and angular displacement 0 . Both w and 0 are 
assumed to be large enough to be significant but small when compared 
with the radius of the ring. Therefore, all terms can be neglected 
containing second or higher powers or products of w/r or- $4 Also 
all terms containing second or higher powers or products of the 

following items are neglected: %/r, L M w here r 

Wl d0 2 / & d d0 2 

and 0 are the polar coordinates of the ring. 

This work was conducted at the University of Michigan under the 
sponsorship and with the financial assistance of the National Advisory 
Committee for Aeronautics. 


where 

V p o 

i/Pi 

Mi 

El 


SYMBOIS 

a,b major and minor axis, respectively, of elliptical ring 

a n*' b n> c n> a n' e n> f n'5i' h n Fourier coefficients 



q. applied pressure load per unit length 

r , Q polar coordinates of original ring 

r l , ^l polar coordinates of deflected ring 

r A radius from origin to point A on original ring 


s 


arc length along circumference of ring 
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nondimensional parameter 

nondimensional parameter 

displacements of a point on ring parallel to x- and y-axis 
respectively 

radial deflection of ring 

Fourier coefficient for radial-deflection function 
bending flexibility of ring 
certain functions 


(El/qa 3 ) 
(ei/ qK 3 ) 




height above x-axis of straight-line portion of rings of 
family II 

axial stress at point A on ring 

polar moment of inertia of equivalent ring of elastic weight 
referred to original ring 

polar moment of inertia of equivalent ring of elastic weight 
referred to deflected ring 

certain constants 

bending moment 

bending moment in ring if change of geometric shape of ring 
is neglected 

bending moment at point A on ring 

bending-moment expression derived from consideration of 
applied pressure load and final shape of deflected ring 

bending-moment expression derived from consideration of 
internal stress due to change of curvature of ring 

bending moment in ring at point defined by 0 

certain constants 

radius of circular arc of rings of family II 

area of equivalent ring of elastic weight referred to 
original ring 
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area of equivalent ring of elastic weight referred to 
deflected ring 

a- angle of rotation of cross section of ring 

Pq radius of curvature at a point on original ring 

Pi radius of curvature at point on deflected ring corresponding 

to p Q 

$ angular displacement 


THEORETICAL ANALYSIS 


In the following analysis two expressions for the tending-moment 
distributions are derived. The first expression is obtained from 
consideration of the action of the applied pressure load on the final 
shape of the deflected ring. The second expression is obtained from 
consideration of the internal stress distribution due to change of 
curvature caused by loading of the ring. Hereinafter the bending 
moment corresponding to the first expression is designated by M a 

and the bending moment corresponding to the second expression is 
designated by . By equating the two expressions, a differential- 

integral equation -which represents the equilibrium condition of the 
final deflected ring is obtained- The solution of this differential- 
integral equation determines the radial deflections of the ring and 
the bending-moment distribution. From the radial deflections, the 
angular displacements can be found from the nonextension theory, which 
is treated in detail in appendix A. 


Expression for M a 

The moment M a is the bending moment determined from consideration 

of the applied pressure load and the final deflected position of the 
ring. Before the expression for Mg is derived, a typical method 

which has been used to determine the bending moment in the ring without 
consideration of the change of geometrical shape of the ring is 
discussed- 

Let the shape of the given ring be defined by 


= r(e) 


r 
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The ring is assumed to be symmetrical with respect to both the x- and 
y-axes. This simplification does not affect the validity of the method, 
which can be used for any shape. The intensity of the internal pressure 
on the ring is designated by q. . 

An imaginary cut is assumed at point A (see fig. 2) and two 
unknowns H A and M A are introduced, where H A is the axial force 

at A (positive if in tension) and M A is the bending moment (positive 

if the outside fiber is under compression) • There is no transverse 
shear at point A because of symmetry. (See reference 2.) The ring 
is now statically detennined. 

The bending moment at any point C, defined by 0, as shown in 
figure 2, can be expressed as 



Substituting the expression 



into the foregoing equation, there results 


Mg = A + Brg sin 0 + Crg 


2 


( 1 ) 


where 



B = H A - dr A 


( 2 ) 



2 


The terms A, B, and C of equation (2) are independent of 0. 


NACA TN No. 1692 


7 


2A 

* 

The angle of rotation a of the cross section at C and the horizontal 
displacement u of point C, "both relative to point A, are given by the 
equations: 



( 3 ) 



sin ©) ds 


(4) 


Here C can be any point on the ring. If A' denotes the other end 
of the ring at the cut, equations ( 3 ) and (4) should also hold at 
point A* ; and since there is no rotation nor horizontal displacement 
at the cut, the following relations are true: 


a A ' “ 


Mg ds 
El 


= 0 


(5) 


u. 


■ • / »0 


r^ - rg sin ©) ds 


M0 r 0 

El 


sin 0 ds = 0 


(6) 


The substitution of equation (l) into equations (5) and (6) yields 
the following equations: 


/W 


r a Sin 0 / 

-2 ds + C <b -2- ds = 0 

El , /El 


A 


El 


(7) 
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r q sin 8 ds 
El 


sin^Q ds 

~ H * c 


Tq- 1 sin 8 ds 


El 


= 0 (8) 


Because of the property of double symmetry, 


r fl sin 8 = 0 

9 El 


r fl 3 sin 9^ = 0 
9 El 


(9) 


Therefore 


A = 


r 0 2 ds 
El 


ds 

El 


B = 0 


( 10 ) 


Or, with the following designation, 


J = 


r 8 ds 
El 


S = (b M 


ds 

El 


(ID 


it follows that 


A = 


-ic 


B = 0 


( 12 ) 
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After substituting equation (12) into equation (l), the expression 
for the bending moment becomes 


The expressions for J and S can be easily memorized because they 
are equivalent to the polar moment of inertia and area of the 
corresponding ring of elastic weight, respectively. (See appendix B 
for more details.) 

It should be noticed that the bending moment given by equation (13) 
is a function of r. Any change of shape of the ring changes the values 
of r and consequently the bending moment also is different . All 
papers in the past have neglected this change and have called Mg, 

given by equation (13), the final bending-moment distribution. It is 
shown later in the examples that the error is considerable for flexible 
rings. The derivation of equation (13), however, leads to the 
establishment of the expression for M a , which is based on the final 

deflected shape of the ring rather than the original shape . 

The final position of the deflected ring is defined by r^ = r^(0). 

Because of the assumption that the ring is inexPensible, or in other 
words, ds = Constant, there Immediately follows from equation (13) 
the following expression for M 0 which is referred to the final 

a 

deflected ring: 




(13) 



(14) 


where 



(15) 
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Equation ( 14) gives the expression for the bending-moment distribution 
from consideration of the applied pressure load and referred to the 
final deflected position of the ring. 

If the shape of the final deflected position of the ring is 
known, the bending -moment distribution can be obtained immediately 
from equation (14). The problem now is to find a method for deter- 
mining the final deflected position of the ring, which, in turn, 
depends on the bending -moment distribution. 


Expression for Mj_ 

The moment Mj_ is the bending moment built up from the fiber 

stresses resulting from the changes of curvature of the ring when it 
is deflected. To set up the relationship existing between the bending 
moment M-^ and the deflections is rather difficult, especially when 
deflections in both directions (two-dimensional) are to be considered* 
Timoshenko presents the derivation of a differential equation that 
gives the relationship between the radial deflections and the internal 
bending moment. (See reference 3*) The differential equation is 




where w is the radial deflection of the ring. This equation is based 
on circular rings and can be used only for rings that are nearly 
circular. In the following paragraphs, a differential equation is 
derived which is more general than the one given by Timoshenko. 

Referring to figure 3, let the original shape of the ring be 
defined by r = r(0) . An arbitrary point A, after loading, experiences 
an angular rotation 0 and a radial deflection w = BD. Its new 
position is completely defined by 0 and w, both of which are 
functions of 9 . The deflected shape of the ring is then represented 
by the equation 


r 


1 " 



(16) 


where 


1*1 = r + w 
0 -^ = 04-0 


(17) 
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From the calculus, the original curvature of the ring at point A is 
given by 



(18) 


and the curvature of the deflected ring C-^ at corresponding point D 
Is given by 



Since 


ie (r + v) r 1 + w' 
^1 ~ A(e + (jf) TTJ 7 


( 20 ) 



r" + w" _ r ' + V 1 0 M 

(1 + < f >') 2 U + 0') 3 


( 21 ) 


where r', w', and so forth indicate derivatives with respect to 6, 
equation (19) becomes 
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(r + w) 


2 


(r + v) — — i— £ E l* v' 0" 

L(n-ti') g (1.1*'] 3 


<- + V>2 * ('tT^ 


3/2 



( 22 ) 


From the assumption that the ring is inextensihle, there exist 
certain relationships between the angular displacements 0 and the 
radial deflections w. These relationships, as given in equations (23) 
and (2k), are derived from the nonextension theory. (See appendix A.) 



(23) 


r = 






(2k) 


The foregoing relations, when applied to circular rings for 
which r = Constant and r' = 0, become 


d' = -2 

r 



(25) 


which are the equations given by Timoshenko (see reference 3 , P* 208) 
for inextensihle circular rings. 
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Substituting equations (23) and. (2k) into equation (22) and. 
simplifying (see appendix C for details), the following equation is 
ob tained: 




Now, from equation ( 18 ), 


( 26 ) 



Therefore 


Pi 




(27) 


But as already mentioned (assumption (2) of the INTRODUCTION) the 
following equation is true: 




El 
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Therefore 



(28) 


This equation expresses the relationship between the bending moment 
at a certain point on the ring and the radial deflection of the ring 
at the same point- 

Equation (28) is valid not only for rings but can be applied to 
any curved bars provided that the contour of the bar is regular and 
smooth, the deflections are not too large, and the bar is inextensible . 
Two extreme cases are given in the following paragraphs for 
illustration. 

First, for circular rings, r is constant and r' = 0- 
Equation (28) reduces to 


which is exactly the formula given by Timoshenko. 

Next, consider a straight beam. Let the origin be chosen at 
infinity so that r = 00 and r ' = 0 at all points on the beam. 
Equation (28) can be simplified as 


For the x,y -coordinate system used in beam deflection, w = y, 
r dJ9 = dx. (See fig. 4.) Thus the foregoing equation becomes 


If 

V + V = 




which is the familiar straight-beam formula . 


A 
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Differential-Integral Equation 

Equation (l 4 ) gives the expression for Mg which is determined 

by considering the applied pressure load and the final deflected 
position of the ring. Equation (28) gives the expression for M^. 

It is determined by considering the change of curvature of the ring 
resulting from the loading. At the final deflected position of the 
ring, these two bending moments Mg and should be equal at every 

point on the ring; that is, 

M a = Mi 

Therefore, 

♦ * - i + (t ) 2 


In this equation, r^ and are referred to the final deflected 

position. (See equation (l 4 ).) 

From equation (17), 


r-i = r + w 



(r + w)^ 



(30) 


Therefore 


J 


1 " 





1+22 

r 


+ 



(3D 


Remembering the assumption that any term containing a power 
of w/r higher than second can be neglected, equations (30) and (31) 
become 



+ 2rw 


( 3 ^) 


1 6 
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and 


J i =i r2 S + i ** i 


= J + 4 2wr — 

.1 El 


(33) 


Putting these expressions into equation (29) yields 


v" + w = — — (r^ + 2rw 
2 El V 


-§-§/«• tf )^ 1 + (t/ < 3,1 > 


Transferring all terms containing w to one side of the equation gives 


v 


" + w - + Ct ) 2 * £ + (t ) 2 / ”• g 


1 

2 El 


V 1 + (t ) 2 - HK# 


(35) 


This is the final differential-integral equation with the radial 
deflection w as the only variable in the equation. The solution of 
this differential equation and a method to make it more practical for 
applications is discussed in the following section. 


Solution of Differential-Integral Equation 

The solution of the foregoing differential-integral equation can 
be made by the following steps: 

(a) All the known functions and unknown functions are expanded 
into Fourier series, with known and unknown Fourier coefficients. 


r 
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(b) By comparing the Fourier coefficients of both sides of the 
equation, a system of simultaneous equations is obtained. 

(c) The solution of the simultaneous equations gives the unknown 
Fourier coefficients. 

In carrying out these steps, rewrite equation (35 ) , using the 
following relationship from the elementary calculus: 



Equation (35) then becomes 



In equation (37 ) , let 
functions : 


F-j_, F 2 , and 


stand for the following 
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The functions F^, Fg, and F^ are known for any given ring. The 
equation becomes 


v" + w - wqF 2 + 


F x w d 9 



(39) 


Since F^_, Fg, and F^ are known functions, they can be expanded 
into Fourier series with known coefficients. Therefore, 


F 1 = 2_ a n cos n0 

n=0 ,2,4. . . 


f 2 = 


b n cos n 9 


(40) 


n=0 , 2 , • • • 


Fo = 


c n cos n G 


n=0,2,4. . . 


where a n , b n , and c n are known Fourier coefficients. For double 
symmetry, the case under study herein, n is always an even number. 

Now let w be represented by a Fourier series with unknown 
coefficients : 


w 


Hfc:0,2,4. . . 


A m cos jdB 


(41) 


where A m 's are to be determined. From the fact that the ring is 

double-symmetrical and the external load system is also double- 
symmetrical, the radial deflection w also is a double -symmetrical 
function. Therefore all m's are even numbers- 
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Differentiating equation (4l) twice gives 


v" = - cos m0 

m=0,2,4. . . 


(42) 


Therefore, 


v + v = 


(l - m 2 ) 


A m cos m 6 


m=0,2,4. • • 

The substitution of these relations into equation (37) yields 


(43) 


H (! - m^A^ cos m0 - q ^ b n cos n0 ^ A^cos m 8 

m=0,2,4... n=0,2,4... m=0,2,4. 


+ * XT 

S n=0,2,4. . . 


a n cos m0 


A m cos m 6 £ a n cos n0 d 9 

m=0 , 2 , • • • n=0 , 2, ^4- • • • 


V / c n cos n0 - — 

n=0,2 , 4. . . ® 


n=0,2,4 . . . 


a n cos n0] 


m 


In order to facilitate the comparison of the coefficients on both 
sides, equation (44) is simplified. First, consider terms involving 
the integral sign in equation (44) as follows: 


A m cos m0 /_ cos n0 d0 

mt=0,2,4... n=0,2,4... 


J )2« 

cos m 6 
0 n=0,2,4 . . . 


a n cos nd d0 


= + JZ ^Am a m 

m=2,4. . . 


(45) 
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and equation (MO becomes 


k k»k cos mB - q cos n0 k A m cos a £ 

m=0,2,4... n=0,2,4... m=0,2,4... 


+ a 


n=0 


)> a n cos n9 (2rtAQa 0 + ^ nA^a^ j 

,2,4... \ k=2,4... / 


§' 


^ 11— 0 y 2 , ^ 4 - • • • 


c-n cos n0 - — 
n o 


a n cos n0 


XI — 0 y 2 , ^4" • • • 


(46) 


The second step is to simplify the second term of the equation 
which contains the multiplication of two Fourier series. By actually 
carrying out the multiplication, 


t> n cos nB 


cos m B 


n — 0 , 2 , ^ l- • • • 


nt=0,2,4 . . 


= + d 2 cos 2 6 + d^ cos kG + 


= ]> <V COS me (47) 

m=0 ,2,4. • • 
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where 


<*0 = + \ XI 

* n=0/2~£. . . 

^ = |(Ao b 2 + A2 b o) + | XZ (An b n +2 + An+2 b n) 


Li-0 y 2 y U • • • 


d 4 = f(Ao b 4 + V 2 + A U b o) + § o XI (An b n+4 + An+4 b n) 


(48) 


V - + ^^-2 + • • • + A m _ 2 h 2 + 


+ 


1 

2 


H 

n=0,2,4. . . 


(An b n+m + An+mhn^) 


and equation (46) then becomes 


(l - m*") cos me - q cos me 

3 y ^ • • • m=0 y 2 y ^ • • • 

+ ^l 2 ^ 8 ^ + Z ^A^k) Z % cos me 

\ k=2,4... /m=0,2,4... 


4 

2 


m=0,2,4. . . 


^m ‘ | a m) cos 1110 


(49) 
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or 


1 


ife0,2,4 . 


(l - ^ ( 2«Aq3q + -nA^a^^ 

\ k=2,4... / 


cos m£) 


lx... ('’“"s 8 ”) 


2 m=0,2,4. • • 


cos m0 


(50) 


Equation (44) has been reduced to a more usable form in equation (50). 

Now the Fourier coefficients on both sides of equation (50) can 
be compared. 


For m = 0, 


Ao " 4^ + | ( + 


rtA \k a ln a 0 = §(c 0 “ f®o) 


k=2,4 . . . / 


For m = 2, 


-3^2 - qdg + ^ (2^80 


♦ I 


k=2,4. 


’'A-k a k) a 2 


= f (°2 • jh) 


(51) 


For m = 4, 


- qd^ + | ^nAQao + = ^(04 - 

\ k=2 , 4 . . . / 


As many simultaneous equations as wanted can be formed. The - solution 
of m simultaneous equations gives m unknowns, Aq, Ag, A^ . . . . 

In actual cases, as shown in the examples in appendix E, three or 
four simultaneous equations, which give three or four A m terms, are 

sufficiently accurate for ordinary rings • The following system is the 
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system of four simultaneous equations -which give four Fourier 
coefficients, Aq, A 4 , and Ag. The following equations are 

obtained from equations (5l) and (48) after collecting the terms: 


,(l - qb 0 + 2§ 


+ 2£rta 0 a 0 ) + A2(|b 2 + *§a 0 a 2 ) + A 4 (-§b 4 + rt|a 0 a 4 ) 


+ A 6(-| b 6 + §* a 0 a 6) = |(°0 ‘ g a o) 

Ao(-l b 2 + 2,1 f a 0 a 2) + ^(“3 ' l b 0 ' f b 4 + 7t g a 2 a 2^ 

+ A4^|bg - |b 6 + + A 6(-| b 4 - f b 8 + I | a 2 a 6) = |( c 2 ' 

Ao(-qb 4 + S^aQa^ + Ag^bg - |bg + flga^a^ 

+ A 4 (-15 - qb 0 - |b 8 + r|a 4 a 4 ) + Ag^bg - |b 1Q + *§a 4 a 6 ) 


K52) 


= f ( C!+ " I 84 ) 


Ao^-qbg + 2 ff |a 0 a 6 ^ + Ag^|b 4 - |b Q + i^agag) 

+ A 4(' 2 b 2 - § b 10 + I f a 6 a 4) + ^("35 “ l b 0 ' 2 b 12 + T f a 6 a 6) 
= |( C6 ' s a 6 ) 
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With A m determined., the radial deflections can he obtained from 

equation (4l) and the bending-moment distributions, from equation (14) . 
The angular displacements 0 can be determined by equations from the 
nonextension theory (appendix A) . 


Angular Displacements 

When the radial deflections w are known, the angular displace- 
ments can be found from the following equation, obtained from the 
nonextension theory (appendix A) • 


r = 


_ w _ 

r r 


ll 

r 


(53) 


The integration of equation (53) gives the angular displacements 0. 


0 = 




(54) 


Bef ore 
cients e n 


evaluating the foregoing integration, the Fourier coeffi- 
for the function l/r must be determined: 


1 

r 


z 

n=0,2,4 . . . 


e n cos n0 


(55) 


Since l/r is also a double-symmetrical function, the n’s are even 
numbers only. Therefore, 


w 

r 


1 


afcO ,2,4... 


A m cos m 6 


I 


e n cos n 9 


ib=0 . 2 . • • • 


- z 


f n cos n 0 


( 56 ) 
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where 




= k*0°0 + i XI Vfc 

n=0,2,4... 


f 2 = o( A O e 2 + A 2 e o) + o y ( A n e n+2 + A n+2 e n) 

1 ^ 0 , 2 , 4 ... x y 


> 


(57) 


= ^(Ao e 4 + ^2 e 2 + A 4 e o) + k y ~ ( A n e n+4 + A n+U e n) 
2V n=0,2,4... 


and 


2 de = 


(f“o + fp 008 + ^4 008 ^ + • • 


d0 


f 0 e + 


n ^2 ^ ^ 4 - * • 


— sin n 0 
n 


(58) 


It is also known that 


t - 

r n=2,4... 


sin nd 


(59) 


and 


r r 


mAj^ sin m 9 


ne n sin n 9 


nt=0,2,4. . . 


n=2,4 . 


rt= 0 , 2 , 4 . . . 


g n cos n9 


(60) 
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where 


V - 2 

Sb = g 2_ n A n e n 
ru=2 , 4 . . . 


82 = 


n=2,k . . . 


n(n + 2)(A n e n+2 + A n+2 e n ) 


( 61 ) 


84 = \ 22 n ( n + 4 )( A n e n+4 + Vt-4 e n) “ (&%) 

11=2 y ^ 4 - • • • 


Therefore, 



cos n0 dj9 


= 'Bo 9 ' 22 

n=2,4. . . 


8n 

n 


sin i)9 


( 62 ) 


Substituting equations ( 58 ) and (62) into equation (54), the 
following form for (f is obtained: 


0 = ' f n» 


sin n0 - g^Q - 


n=2,4 . . . 


n=2 , • • • 


— sin n0 ] + C 
n 7 


fe-gb) e + 21 — 


- «n . 


n=2 y ^4- • • • 


n 


sin n Q 


+ C 


( 63 ) 


where C is the constant of integration. 
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The present boundary conditions require that: 


(1) at 0=0, = 0 

and (2) at 0=2 0 = 0 


( 64 ) 


For condition (l), equation (63) becomes 


0 ( 0 ) =0 = 0+ 2_ (°) + C 

n=2 , 4 . . . 


C = 0 


For condition (2), equation (63) becomes 


0(f) « 0 - (f 0 - 8b)f - 51 (0) 

W 2 n=2,4... ■ 


( f 0 " ^0)2 ” ° 
f o - fib = 0 


And equation (63) becomes 


0 = 


n=2,4. 


fn ' 811 3 in n8 
n 


m 


which is the final form for the angular displacement 0. 
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The fact that fQ - gQ must he zero gives a good check on the 
values of Aq, Aq, Ag . • • • This condition should always he 
satisfied. A proof that this condition fQ - ^ = 0 is always 
satisfied automatically is given as follows: From equation (5), 



HI* 

p 

OQ 

II 

O 


whereas from equation (28) 


M w + w 



V \ r / 


Therefore, 


i 

M rts l w + w " ds 

EI M-e/ 

Pin 

_ / w + w d0 = 0 (66) 

Uo r 


Since 


1 

= y e n cos n9 

r 

rtO ,2,4... 

ft 

w + w = 

\ (l - m 2 ) A m cos me 


m=0 ,2, 4. • . 


therefore 


w + w" y 

r m=0,2,4 . . • 

fl - rn^Arri cos me zL e n cos n9 

V 7 n=0 ,2,4... 

- H 

n=0 ,2,4. . . 

h n cos n.6 


cos n 0 
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viiere 


*0 




hg - 


(67) 


Substituting into equation (66) gives 



= 2ith Q = 0 


or 


ho = 0 


Comparing equations (57)> (6l), and ( 67 ), it can be seen that 


^0 = f 0 " &0 = 0 


Thus fQ - 8 q = 0 is a condition which must be satisfied automatically . 


Summary of Procedure 

In the foregoing sections, there has been developed an analytic 
method of finding the bending-moment distribution in double -symmetrical 
rings of arbitrary shapes, acted upon by internal pressure loads. The 
ring must have a regular and smooth contour, with its thickness 
negligible compared to the radius . For such rings the procedures to 
solve this problem can be summarized in the following steps: 
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(l) Determine the Fourier coefficients of the following given 
functions: 


' x y n=0,2,4. . . 


cos n/9 


F 2 = 


El 




n=0,2,4. . . 


b n cos n 9 


f, = d\ 1 + fcf = 


3 21"^ = n=0,£t... * 


c„ cos n0 


( 68 ) 


If the contour of the ring is expressed analytically in simple functions, 
the Fourier coefficients can he obtained by the usual method of 
integration. Otherwise, the coefficients must be obtained by graphical 
integration, which is explained in appendix D* 

(2) Perform the following integrations: 


J = 


El 


ds = 




S = 9 — = 


ds 

El 


^2 it 


sf : W ie 


(3) Calculate the following constants by substituting the results 
obtained from steps (l) and (2): 

*00 = “ b 0 + l^O^ 


*02 = 
*04 = 
*06 = 


"kp 


+ a . 

2 + ?0*k 

2 * fw 


A 


NACA TN No. 1692 


31 


Kgo = -b 2 + ^ a 2®0 

*22 = "' b 0 - -jf + f ® 2 ® 2 

bo b£ „ 

*24 = “ " 2 " + s®2 a 4 

*26 = ’1" " ~2 + s® 2 ® 6 

*40 = _b 4 + ^ Ia 4 a O 
*42 = ~~§ " Y + | a 4 a 2 

K 44 = ^0 - T + S® W 

xr _ b 2 _ b 10 . *. - 

k 46 2 ” “ 2 “ s^6 

k 6o = " b 6 + ^6^ 

"bi bo . 

*& ' 2 2 + I® 602 

%*=-!• ^*1*6*4 
e 66 = - b 0 * ¥ + s a 6 a 6 

«o - K° o i"o) 

«2 - f (=2 - l^) 

41* = f(°4 - |*1.) 

«6 = 5 (=6 - f*6) (69) 
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(*0 Substitute the foregoing constant terms into the following 
system of simultaneous equations: 


A cQ l + *(X>) + A 2 K 02 + a 4 K 04 + A 3k>6 = % 
A^Kgo + A 2^q + ^22^ + A 4^24 + -^ 6^26 = ^2 
Ao k 4o + a 2 k 42 + a 4^^ + K 4t) + a 6^46 = ^4 
a O K 6o + ^62 + a 4 k 64 + A 6^^ + ^66) = % 


(70) 


(5) For a given pressure load q, the foregoing simultaneous 
equations can be solved for Aq, A^, A^, and Ag. 

(6) The radial deflections w are given by the following equation: 


w = Aq + A^ cos 20 + cos 40 + Ag cos 60 (71) 


(7) In order to find the angular displacements, the Fourier 
coefficient e n of the following functions should be first determined: 


7 = 2Z e n 008 110 (72) 

n=0 ,2,4. • . 

Then values of A n and e n are substituted into the following 
equation to see if the condition f Q - = 0 is satisfied. 

f 0 “ So = + 2 C 1 " n2 ) A n e n 

n=2,4 . . . 


( 73 ) 
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(8) Knowing e n and A n , determine f n and as given in the 

following equation: 

f 2 = ^( A 0 e 2 + ^ 0 ) + \ 5 __ ( A n e n +2 + A n+ 2 e n) 

f k = i(Ao 0 4 + A2 e 2 + A 4 e o) + 2 5_ (An®n+4 + An+4®n) 

n=0,2,4 . . . 

f 6 = ^( A 0 e 6 + A2 e 4 + A 4 e 2 + Vo) 

+ i 5Z ( A n e n +6 + An+ 6 e n) 

^ d^O.2,4*** 

^ (74) 

82 = i Y n(n + 2 ) fVn +2 + A n+ 2 e n) 

2 n=2,4... 7 

84 = \ Y n ( n + 4) (Vn+4 + A n + 4 e n) ' ^2*2 
2 n=2,4... ' 

86 = | 51 n < n + 6 )(Vn+6 + V6 e n) ‘ ^ ( ft 2 e 4 + V&) 

n=2 , 4 . • • 


( 9 ) The angular displacements 0 can be obtained from the 
following equation: 


f-- H 

n=2 y • • • 



n 


sin D0 


( 75 ) 
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(10) The moment distribution M is given by the following 
equation : 






£ + 2rw 
S 




(76) 


From the foregoing procedures it can be seen that the main part of 
the solution of the problem is the determination of the Fourier coeffi- 
cients of several known functions • The rest of the procedures are simply 
algebraic. On the average, it takes 20 hours to solve completely a 
problem. 


Preparation of Charts 

Following the foregoing procedures, solutions have been obtained 
for two specific families of rings. The first family consists of 
elliptical rings (fig. 5(a)) with various eccentricities. The second 
family consists of rings formed by two semicircles and two straight 
lines. (See fig. 5(b).) The ratio of the radius of the semicircle 
to the height of the straight-line portion is variable. All rings are 
assumed to be of constant El. 

Charts and tables are provided for both families of rings. The 
Fourier coefficients a n , bjj, c n , and e n are given in tables 1 and 2 

and are also plotted in figures 6 to 13* Since both families of 
rings are of constant El, the Fourier coefficients are obtained in 
terms of El and the dimensions of the rings. 

Values of J and S are given in table 3 and also plotted 
in figures lk and 15* When these values are substituted into 
equation (70) the simultaneous equations are obtained. For a given 
value of q and El, the simultaneous equations can be solved for A q, 

A 2 > a 6 * ' * * Figures 1 6 to 21 give the values of A n for 
various ratios of q/EI. 
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Knowing A n , w, <f), and M are obtained from equations (71)* (75)* 
and ( 76 ), respectively. The results are plotted in figures 22 to 42. 

Two examples are given in appendix E to illustrate the foregoing 
procedures in detail. 


TEST 

Test Specimens 


Two steel rings were tested, one elliptical and the other made 
up of two semicircles joined by two straight lines. (See fig. 5«) 


Test Apparatus and Procedures 

Figures 43 and 44 show the test setup. The ring was tested in a 
horizontal plane. Its circumference was divided into a number of 
segments with equal arc lengths. Wires, loaded equally, were attached 
to these points. Each wire was led through a pulley, connected to a 
fixed horizontal ring in such a manner that the pull on the test ring 
was normally outward at each point. (See fig. 45-) The concentrated 
loads were sufficiently close so that they could be assumed as 
simulating pressure . 

In order to register the deflections of the ring, a wooden board 
was placed on top of the ring. The contour of the ring was traced on 
the board, with marks denoting the division points, before and after 
the ring was loaded, as shown in figure 46. The radial and angular 
displacements can be measured directly on the board. 

Bending moments were found by electrical strain -gage readings at 
three points on the ring. (See fig. 47-) The value of the effective El 
used was determined from bending tests on a specimen cut from the ring. 

It might be mentioned that the test was at first tried with the 
aid of a pressure bag. The ring was laid around the bag which was 
blown up with compressed air. Special devices were used to adjust the 
pressure between the bag and ring so that the ring would be uniformly 
loaded- The test was unsuccessful, however, since it was found that 
uniform loading could be obtained only through long and tedious work. 

(See fig. 48.) 
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Test Results and Discussion 


The test results are plotted in figures 49 to together with 
curves obtained by the method developed in this paper. The agreement 
between test data and calculated results is good. 

In figures 51 and where bending-moment distributions along the 
circumferences of the rings are given, an additional curve is shown in 
each figure. These added curves represent the bending -moment distri- 
bution when changes of geometric shapes of the rings caused by loading 
are neglected. The symbol Mq is used to designate the bending moment 

calculated without considering the deflections of the ring. The value 
of Mq can be found from the following equation: 


It can be seen from f igures 51 and 54 that the bending moment 
obtained without considering the change of geometric shape is too 
conservative • At the point of maximum bending moment the difference 
is quite large . 


From an analytical method derived for the determination of bending- 
moment distribution and radial and angular displacements of flexible 
rings of arbitrary shape under internal pressure load, with the change 
of geometric shapes caused by the load being considered, the following 
conclusions can be made: 

1. Results obtained by the present method showed good agreement 
with test data and proved that results obtained when change of geometric 
shape was not considered are inadequate for flexible rings. 

2 . Although only rings with double symmetry under internal pressure 
loads are discussed in the present work, the method can be extended to 
include rings of any shape under any system of external loads. 

University of Michigan 

Ann Arbor, Mich., October 17, 19^6 
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APEETOIX A 
NONE.XTENSION TSEOEY 


The deflections of a loaded ring can be completely defined by the 
radial deflections w and angular displacements 0, both of which are 
functions of 6 . (See fig. 3*) The deflections w and (f are 
independent of each other if there are no additional conditions imposed 
on the ring. They are definitely related, however, if the ring is 
assumed to be inextensible • The nonextension theory gives the relation 
between w and 0 for such rings . 

Timoshenko (reference 3> p- 208) has given the relation between 
the radial deflections w and the tangential deflections v (which 
correspond to angular displacements 0 in the present case) for 
circular rings as follows: 


g + - - ° (ai > 

For circular rings this can be written as 

rf> ' + w = 0 (A2) 


where 


d> ' = M 

d.e 


This relation, however, is not quite accurate for rings of noncircular 
shapes • A relationship between w and 0 which is more general and 
can be applied to rings of any shape is derived by means of several 
steps as follows. 


Arc Length in Polar Coordinates 
Given a 3mall element AB = ds on an arc C . 


OA = r 


OB = r + dr 
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With 0 aa the center, OA as the radius, swing a circular arc which 
will cut OB at D. Then 


OD = OA -= r 


BD = dr 


and 


AB = (ds)^ = \/(AD) 2 + (BD) 2 


= y(r d 0)2 + (dr) 2 



where the prime means derivative with respect to 0 • 


(A3) 


Increase of Arc Length Due to Angular Displacements 

When point A (fig. 55) is allowed an angular displacement <f> 
and point B an angular displacement 0 + d^, with no radial deflections, 
the new position of the element AB Becomes EF . Since there is no 
change of radius. 


OE = OA = r 


OF = OB = r + dr 


and the length of the arc EF is 

EF = (ds)^ = \/r 2 (d0 + d0) 2 + (dr) 2 


= r d0\|7l + 0') 2 + (^J 
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The increase of arc length due to the angular displacement is 
where 


(Ads ) 0 


= (ds)gp - (ds)^jg 


= r d0 





(A5) 


Increase of Arc Length Due to Eadial Deflections 

The element of arc EF is now allowed to have radial deflections * 
Point E moves radially to G and F, to H. (See fig* 5 6 .) 


GE = w 


FH = w + dw 


Tne arc length of GH is then 


GH = (ds) GH = \f( GI) 2 + (IH) 2 


= y( w + r)(d0 + d 0) 2 + (dw + dr) 2 


= (r 4- v)deW(i + <t') 2 + (A6) 


and the increase of length of the arc due to radial deflections alone 
is (Ads) w where 


(Ads) v = (r 4- w)d0 \| ( 1 4- fi'Y + 


+ w 


r + w 


-r 


d© V ( 1 + 0') + 



(A7) 
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Equation for Deflections of an Inextensible Eing 

The nonextension theory requires that the total increase of the 
arc length due to both radial deflections and angular displacements 
should be zero. In other words. 


(Ais)0 + (Ads) w = 0 


From equations (A 5 ) and (A 7 ) , 


(r + w)de\ /(l + $'Y 


+ w 


r -1- w 



or 


<r ♦ - rf 7(2 




(a8) 


Squaring both sides and dividing by 


yields 


or 



or 





> 2 + 
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or 


r 



- 1 


(A9) 


Equation (A9) gives the relation between the angular displacements and 
the radial deflections for an inextensible ring- Equation (A9) is 
very general and can be used for rings of any shape- 


Simplif i cations 

Equation (A9) can be greatly simplified for rings, the 
deflections w and 0 of which are not very large when compared 

with the radius of the ring. In other words, ^ « 1 and 0 « 1. 

Also, it must be assumed that ~ « 1 and 0' « 1- Then any term 

of power higher than two or products of the foregoing items can be 
neglected. Equation (A9) becomes 

(l + 2 )( 1 + 20 ') = 1 “ 2 ^.^ 
or 

1 + 2d' + 2^ = 1 - 2*1 vl 

r r r 


or 


r = 



(A10) 


Differentiate equation (AlO) once as follows: 


0 


w r_l 
r r 





(All) 


Equations (AlO) and (All) are the relationship between 0 
inextensible rings. 


and v for 
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For circular rings for which r = Constant 
equations (A10) and (All) reduce to 


r = 


w 

r 


r 


v f 
r 


and r T = 0 , 


which are exactly the relations given hy Timoshenko in reference 3 . 
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APPENDIX B 

EQUIVALENT RINGS OF ELASTIC WEIGHT 


The term "elastic weight" has been used, by Mohr (reference 4). 
It is designated by "dW" and is defined as 


dW = — 
El 


(Bl) 


The equivalent ring of elastic weight is then defined as the ring, 
the median of which has a shape which is exactly the same as the shape 
of the original ring, but the thickness t^_ of the equivalent ring 

is l/EI, ■where El is the bending flexibility of the original ring. 
The cross-sectional area of the equivalent ring then is 

48 - / g <B2) 

and the polar moment of inertia of the equivalent ring is 


I = j/ 7 r 2 t x ds = 




When equations (B2) and (B3) are compared with equation (ll), it is 
seen that 


S = A 
J = I 


of the equivalent ring of elastic weight. 
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APPENDIX C 

DETAILED DERIVATION OF EQUATION (27) 


The details of deriving equation (27) are as follows. Start with 
equation (22), where 


(r + w) - (r + w) 


r" + w' 


r 1 + w 1 


(i + 0') 2 d + 0') 3 


r 


+ 2 


. .\ 2 

r + w ' 


W 7 


(r + + ( T -rrf 2 


3/2 


d?) 2 -( l+ ?) 


r 

r 


M II 


W 

r 


r f v 
r r 


(1 + f ) 2 (1 + /’) 3 




+ 21 


f Y + Wj 

r r 

1 + jz5 


•(1 * *) ; 


/r + x 


,\2 


n3/2 


1 + 


C 1 + 9 


2 \ 1 + 0 'y 


(ci) 


By neglecting all terms containing second or higher power or 
products of the deflection items and using the binomial theorem, 
the following relations are obtained: 



(1 + 0’)2 



+ - 2 & E- 

r r 


rj. + w_l 

-E Z-4" 

(1 + /') 3 



II 





2 -(— 


r \r 




(C 2 ) 


Substitution of equation (C 2 ) into equation (Cl) yields 



VJI 
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The foregoing equation can be further simplified by using 
equations (A10) and (All). The numerator of equation (C2) becomes 



and the denominator becomes 


Denominator = 






r 




Therefore 


_ 1 _ 

pi 



Numerator 



(C4) 
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V7 


or 



But from equation (l 8 ) 



Therefore 


or 



or 
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APPENDIX D 


GRAPHICAL METHOD TO DETERMINE THE FOUPIER 
COEFFICIENTS a^, b^, AND e^ 


For rings, the contours of which cannot be expressed analytically 
as simple functions, the following procedure is a graphical method for 
quick evaluation of the Fourier coefficients a n , b n , c n , and e n of 

the following functions: 


F 


1 = 



n=0 ,2,4. • . 


a^ cos n0 


f 2 

f 3 



XI 

n=0,2,4 . . . 

Yl 


n=0 ,2,4. . . 


b n cos nd 


c n cos n0 


1 

r 


n 

n=0 ,2,4. • • 


e n cos n0 


Table 4 is filled in first. Then row3 7 to 18 are plotted 
against the arc length s, and rows 19 to 22 are plotted against 6 . 
Let E rj y Eg, . . . E22 represent the area under the corresponding 

curves. The Fourier coefficients are 


a 0 = ^7 

of£ 

II 

i 5 

°o * ^15 

e o = ^19 

CO 

11 

CJ 

cs? 

b 2 = ^ 

03 . 4 e i6 

e 2 = ^20 

=*> * ^9 

b„ = iE 13 

!h 

n 

0 

e 4 = ^21 

9 e ■ S 

b 6 - S 

0 

On 

II 

00 

e 6 = ^22 
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APPENDIX E 
EXAMPLES 


Two examples are given in this appendix. A method of checking 
the final results is also presented at the end of this appendix. 


Example 1 

Given an elliptical ring with the following data (see fig. 5(a)) : 
a = 30, b = 25.ll, El = 81,000, and q = 10. 

With the coordinate system shown in figure 5(a), the equation of 
the ellipse can be expressed either by 

r = b \j 1 - k^ sin 2 0 (El) 

or by 

r = a\Jl - k 2 cos 2 p (E2) 

where 

k 2 = 1 - (if (E3) 

\ 8 / 


The general procedures given under Summary of Procedure are 
followed . 


(l) In order to determine the Fourier coefficients v v 
and c n of the following functions 
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the following formulas can be used.: 



where 


c n = 



cos n0 d0 


8 0n = 1 if n = 0 


6 0n = 0 


If n / 0 


Performing the Integration gives 


*0 = 0 . 84334 (^) 

"2= -°- 15059 (i) 

H - 0.00689^J 

b 0 = O.7760o(||j 

b 2 - 0 .20669(||) 

h = o- 0 * 0 ^) 
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(2) In order to evaluate J and S the following integration 
should be performed: 



Substituting in these equations gives 

J = h -87580 (fjj 

3 * 5 ' 78ll,7 (4) 

Therefore 

| = 0.84335a 2 
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(3) When the values of a n , b n , c n , J, and S are substituted 
into equation (69), the following constants are obtained: 


%> - 

k 02 - 0. 03i.33(^) 


*20 - °- o686 \|ij 

*22 " -o-'n 220 ^) 

% = °- 102 ^) 


K^o = -0.01067 (|j 


K ks = 0.10278(||) 
= 0.1033e(|!j 


Qo 


0 .00280 



02 = -O.0630o(|!) 

\ * °- 011l8 (s) 
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(4) When these constants are substituted into equation (70), the 
following equations are obtained: 

Ao(t " 0.00306) + ^(0.03^33) + A^( -0.00536) = 0.00280a 
Ao(0. O 6867 ) + Ag(-3t - 0.77221) + A^(0. 10278) = -0.06300a 
Aq( -0.01067) + A 2 (0. 10278) + A^(-15t + 0.10338) = 0. 01118a (El*-) 

where 



The parameter t is nondimensional and is a measure of the flexibility 
of the ring. In the present example t = 0*3* 

(5) Solving the preceding simultaneous equations gives 


Aq = 0.00503a 
a 2 = 0 . 03778 a 

A^ = -0 . 00167 a 

( 6 ) The radial deflection w is then given by the following 
equation: 

w = Aq + A 2 cos 2 8 + A^ cos 48 

= (0.00503 + 0.03778 cos 20 - O.OOI 67 cos 48) a (E5) 

The values of w/a for various values of the angle 0 are given in 
table 5 and also plotted in figure 24. 
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(7) The Fourier coefficient e n can be obtained from the 
following formula: 


>*/2 


e„ = 


n "C 1 * 6 0n) Jo r 


— cos n0 d0 


Performing the integration gives 


e 0 = 


1-09979 

a 


e 0 = 


0 .09752 
a 


©I. = 


-0 .00217 
a 


Substituting A n and e n into equation (73) gives 

fo — ~ AqSq — 1 * 5Ag©2 — 7 = "0 *00001 = 0 

Therefore the results are checked. 

(8) The coefficients f n and are obtained from equation (74) 

as follows: 

f 2 = 0.04188 
f^ = -0.00001 

g^ = —0 .00100 

g k = -0.00737 


NACA TN No. 1692 


55 


(9) The angular displacement <f> is, therefore, 

0 = '|( f 2 " 6 e ) sin 20 " ij ;( f 4 ' 84 ) sin ^ 0 

= -0 .02144 sin 26-0 .00184 sin 46 (E6) 


Values of <f> are given in table 6 and are also plotted in figure 31* 

(10) Finally equation (7 6) is used to calculate the bending- 
moment distribution. 


Since 


M = — 


r 2 - — + 2rw 

S 


2 

S 





0 .00302a 2 


and the bending-moment distribution if the change of geometric shape 
of ring is neglected is 


*0 



the bending-moment distribution M can be expressed as 

2 


M 


Mq + qrw - 0.00151qa 


or 


M 


% 


qa qa 


(lh . 0.00151) 

\a a / 


(E7) 


Values of M q /q.a 2 and M /qa 2 are given in table 7 ) M /qa 2 is also 

plotted in figure 38. The last term of equation (E7) is evidently the 
correction needed when the change of geometric shape is considered. In 
this particular case the correction term amounts to about 37 percent at 
the point of maximum bending moment. 
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Example 2 

Given a ring as shown in figure 5(b), with the following data: 

R = 15, H = 9, q = 10, and El = 8k, 375- 

With the coordinate system as shown in figure 5(h), the equations 
for the ring outline are 


where 


^ = k' cos 0 + yl - k'^ sin^0 

r = 1 

R sin 0 


0 <C 0^ a 



(e8) 


k' 


H 

R 


The procedure as followed for example 1 is also used for example 2 
by means of the following steps: 

(l) The Fourier coefficients are 


ao = 1-85526^) 

= 0.78339(|D 
ai, = -0.09307(||) 

b o - 2 -55170(1!) 
b 2 - 1.57561(1!) 
H = -0.02956^ 
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C 0 . 3 . 585^2! Bj 


o 2 = 2.8314-rilj 


=4 - °- 159 Mii, 


(2) The values of J and S are 


J = 

S = 

J = 
S 



16.03321 

8 - 6832 (lf) 

1.8464R 2 


(3) The constants are 


*00 = 

-0 . 06107) 


Kq 2 = 

-0 .26192) 

w 

*04 = 

-0 .04769) 

'r 3 \ 

w 

K 20 = 

-0.52399) 

w 

K 22 = 

-2.31488) 

» 


= -0 .814171 


- -0-°9538(|T) 


k 40 

K kk = -2.54857 


-0.8l4l7(g 
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Qo = 0.07993(|jj 
^. 0 . 69251 ^) 

% - o- 1 ® 6 ^) 


(4) Introduce the nondlmensional parameter t ' = which is 

equal to 2*5 for the present example. The following simultaneous 
equations are obtained: 


2-^3893Ao - 0 . 26196 A 2 - 0 . 04769 A 4 = O. 07993 R 

-0.52394AQ - 9 . 81488 A 2 - 0.81417A4 = 0.69251R 

- 0 . 09538 AQ - 0 . 8 l 4 l 7 A 2 - 4o.o4857A^ = O. 16565 R 

(5) The solution of the foregoing equations is 

Aq = 0.02503R 
Ag = -0 .O7166B 
A 4 = -0.00274E 

(6) The radial deflection w is 

w = (0.02503 - 0.07166 cos 20 - 0.00274 cos 40 )R (E9) 


Values of w/R are given in table 8 and also plotted against 0 in 
figure 28. 
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(7) The Fourier coefficients &Q are obtained as 

e n = Q -7.8538 

0 E 

e = rO - 1 882 9 
2 E 

e 4 = 0-Qi°7Q 

In order to check the results, equation (73) is used. 

f o - gQ = -0.00005 = 0 

(8) The coefficients f n and are then evaluated. 

f 2 = -0.06173 
h f 4 = 0.00537 

82 = -0.00672 

= - 0.02698 


(9) The angular displacements can be expressed by the following 
equation: 

= 0.02751 sin 2 9 - O.OO809 sin k 6 (E9) 


Values of <f) are given in table 9 and also plotted in figure 35. 

(10) The bending-moment distribution is then given by the 
following equation: 



*0 + r w 

qE 2 R R 


0 .01339 


(E10) 


6o 
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Values of M./qRp and Mq jqR^ are given in table 10 and also plotted 

in figure 57 for comparison. The maximum error introduced by 
neglecting the change of geometric shape of the ring is about 35 percent. 


Method of Checking 

A method of checking the final bending moment and displacements 
is given. Example 2 is used for illustration. 

Assume the results obtained from example 2 to be correct. Then 


v = (0 .02503 - 0 .07166 cos 


$ = 0.02751 sin 2 6 - 


M 


Mq 


qR qR c 


+ r w 

R R 


20 - 0.00274 cos 40 )R 
0.00809 sin 40 
- 0.01339 


Using the bending-moment expression, the deflections u and v, 
parallel to the x- and y-axis, respectively, can be determined from 
the original shape of the ring. The deflected ring defined by the 
displacements u and v should agree with the deflected ring defined 
by w and 0 • 


The equations for the determination of u and v are 


Pe no 

Mf, I M„ 

vu)= — r n cos a dfl + ro cos 0 / — = ds 

^ JO El a y Jo El 


(Ell) 


)0 no 

v a = / sin a ds - r^ sin 0 | ^ ds 


'0 


*/2 


El 


a 


!*/2 


El 


(E12) 


Values of u and v, obtained from equations (Ell) and (E12), 
are given in table 11 and are also plotted in figure 5 8 , together with 
the deflections defined by w and 0. The agreement is very good. 
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TABLE 1 .- FOURIER COEFFICIENTS a^, t> n , Cn, AND ^ FOR 
VARIOUS VALUES OF PARAMETER k 

jFamily I, elliptical ringsTj 


OJ. ^ 

1 

H 

II 

0 

0.1 

0.2 

0.3 

0.4 


1 .00000 

0. 94934 

0 .89722 

0 .84334 

0 .78735 

ag^El/a^) 

0 

-0 .05002 

-0.10016 

- 0.15059 

-0 .20160 

a^El/a*) 

0 

0 .00066 

0 .00282 

0 .00689 

0.01354 

/ rnT N 

^2. 

H 

O 

r£> 

1 .00000 

0 .92514 

0 .85052 

0 .77601 

0 .70146 

b 2 (El/a 3 ) 

0 

-0 .07308 

- 0.14210 

-0 .20669 

-0 .26638 

\{sila3) 

0 

0 .00177 

0 .00731 

0 .01703 

0.03151 

o° 

/ t?T s 

H 

"aT** 

1 .00000 

0.90187 

0 .80749 

0 .71684 

0 .62988 

^(ei/s 4 ) 

0 

-0 .09493 

-0.17944 

-0 .25300 

-0.31494 


0 

0.00313 

0.01251 

0 .02817 

0.05014 

© Q ( a ) 

1.00000 

1.02722 

1.05984 

1.09979 

1.15013 

e 2 (a) 

0 

0 .02705 

0 .05899 

0 .09752 

0.14519 

ej+Ca) 

0 

-0.00018 

-0 .00082 


-0 .00217 

-0 .00470 
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TABLE 2 -- FOURIER COEFFICIENTS a^ b n , c n , AND e n FOE 
VARIOUS VALDES OF PARAMETER k' 


(Family II, rings formed, by two semicircles 
and two straight linesT] 


k' = - 
R 

0 

0.2 

o.l 

0.6 

*0 W* 2 ) 

1.00000 

1.261-56 

1.51962 

1.85526 

a 2 (si/R2) 

0 

0.19985 

0.16057 

0.78339 

H( EI / R2 ) 

0 

-0.03398 

-0 .06703 

-0 .09307 

O' 

0 

H 

1.00000 

1.12356 

1.93751 

2.55170 


0 

0.33803 

0.85619 

1-57561 

hW» 3 ) 

0 

-0 .01195 

-0 .06725 

-0 .02956 

O 

O 

H 

15 " 

1.00000 

1.60789 

2.15013 

3-58512 

H 

CVJ 

0 

0 

0.50563 

1.11205 

2.83117 

c^(ei/r 4 ) 

0 

-0 .05622 

-0 .03567 

0.15917 

e 0 (E) 

1.00000 

0.89756 

0.83083 

0.78538 

e 2 (R) 

0 

-0 .07617 

-0.13812 

-0.18829 

e^(E) 

0 

0.01605 

0 .02677 

0 .03070 



TABLE 3.- VALDES OF J AND S 


£ 


k 2 = 1 






Family I 


0.1 


0.2 


0-3 


0.4 


*\ 


j(ei /a-7 
S(El/a) 


6 .2832 
6 .2832 
1.0000 


5.81282 

6.12302 

0.94933 


5-34397 

5.95614 

0.89722 


4.87580 

5.78147 

0.84335 


4.40741 

5-59762 

0.78737 


Family II 


k* - s 

E 

0 

0.2 

0.4 

0.6 

j(ei/e 3 ) 

6.2832 

8.9452 

12.1736 

16 .0332 

s(ei/r) 

6.2832 

7 .0832 

7.8832 

8.6832 


1 .0000 

1.2628 

1.5442 

1.8464 

S \B 7 
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TABLE 4 . - GRAPHICAL METHOD TO DETERMINE THE FOURIER 
COEFFICIENTS a n , Cq, AND ^ 


[see appendix D-J 


Rov 

Item 

Procedure 


Results 

at - 

1 

0 


0° 

10° 

20° 

30° 

40° 

50° 

60° 

70° 

80° 

90° 

2 

r 












3 

m 












4 

cob 2 0 












5 

COB ^ 












6 

COS 60 












7 

r/EI 

(2)/(3) 











8 

(r/El) cob 20 

(7) x (4) 











9 

(r/EI) cob he 

(7) x (5) 











ID 

(r/EI) cob 60 

(7) x (6) 











11 

i^/eI 

(2) x (7) 











12 

(r^/El) cos 20 

(11) x (4) 











13 

(r^/EI) cob he 

( 11 ) x (5) 











14 

(r^/El) cob 60 

(11) x (6) 











15 

r 3 /EI 

(2) x (11) 











16 

(r 3 /El) cos 20 

(15) x (4) 











17 

(r^/El) cos he 

(15) x (5) 











18 

(r^/EI) cob 60 

(15) x (6) 











19 

l/ r 

1/(2) 











20 

(l/r) cos 20 

(19) x (4) 











21 

(l/r) cos 40 

(19) x (5) 











22 

(l/r) cos 60 

(19) x (6) 











23 

Arc length, S 














TABLE 5-- RADIAL DEFLECTIONS v/a AT VARIOUS VALUES OF 9 


[Family I. k 2 = 0-3; t = 0.3J 


CT\ 

On 


1 

2 

3 

4 

5 

6 

7 

9 

(deg) 

cos 29 

cos 49 

0.03778 X (2) 

-O.OOI67 x (3) 

0.00503 

w/a 

(4) + (5) + (6) 

0 

1 .00000 

1 .00000 

0 .03778 

-0 .00167 

0.00503 

0 .04114 

10 

.93969 

.76604 

•03550 

- .00128 

.00503 

.03925 

20 

.76604 

.17365 

.02894 

- .00029 

.00503 

.03368 

30 

.50000 

- .50000 

.01889 

.00084 

.00503 

.02476 

4o 

.17365 

-.93969 

. OO656 

.00157 

.00503 

.01316 

50 

-•17365 

-.93969 

-.OO656 

.00157 

.00503 

.00004 

60 

-.50000 

-.50000 

- .01889 

.00084 

.00503 

-.01302 

70 

- .76604 

.17365 

- .02894 

- .00029 

.00503 

- .02420 

80 

-.93969 

.76604 

- .03550 

-.00128 

.00503 

- .03175 

90 

- 1.00000 

1 .00000 

-.03778 

- .00167 

.00503 

-.03442 
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TABLE 6.- ANGULAR DISPLACEMENTS $ AT VARIOUS VALDES OF 0 


\ 




|Fam1.1y 


I. 


k 2 = 0 -3; 



1 

2 

3 

4 

5 

6 

7 

0 

(4eg) 

sin 20 

sin 40 

-0.02144 x (2) 

-0.00184 x (3) 

f > = (4) + (5) 

( rad . ) 

0 

(deg) 

0 

0 

0 

0 

0 

0 

0 

10 

•34202 

.64279 

-.00733 

- .00118 

- .00851 

-.49 

20 

.64279 

•98481 

- .01378 

- .00181 

-.01559 

-.89 

30 

.86603 

.86603 

- .01857 

- .00159 

-.02016 

-l.l6 

40 

•98481 

•34202 

- .02111 

- .00063 

- .02174 

-1.25 

50 

•98481 

-.34202 

- .02111 

.00063 

- .02048 

-1.17 

60 

.86603 

-.86603 

- .01857 

.00159 

-.OI698 

-•97 

70 

.64279 

-•98481 

- .01378 

.00181 

-.01197 

-.69 

80 

• 34202 

-.64279 

- -00733 

.00118 

- .00615 

-•35 

90 

0 

0 

0 

0 

0 

0 



On 
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TABLE 7 BENDING -MOMENT DISTRIBUTION 


00 


Family I. 


if 


k 2 = 0.3; t = 0.3; § \ = 0.84335- 


1 

2 

3 

4 

5 

6 

7 

8 

0 

(deg) 

r 

a 

(if 

/r\ 2 J 1 

W s a 2 

= (3) - 0.84335 

% 1 
<1 a 2 

= |(M 

V 

a 

from table 5 

v r 
a a 

M 1 
la 2 

= (5) + (7) - 0.00151 

0 

0 .83666 

0 .70000 

-0.14335 

-0 .07168 

0 .04114 

0.03442 

-0 .03877 

10 

.84040 

.70639 

- .13696 

- .06848 

•03925 

.03299 

- .03750 

20 

•85173 

•72545 

-.11790 

- .05895 

.03368 

.02869 

- -03177 

30 

.86989 

•75675 

- .08660 

- .04330 

.02476 

.02154 

- .02330 

40 

.89394 

•79904 

- .04431 

- .02216 

.01316 

.01176 

- .01189 

50 

.92172 

.84956 

.00621 

.00310 

.00004 

.00004 

- .00163 

60 

.95038 

.90322 

•05987 

.02994 

- .01302 

-.01237 

.01606 

70 

•97584 

.95226 

.10891 

.05446 

- .02420 

- .02362 

•02933 

80 

•99365 

.98725 

•14390 

•07195 

- -03175 

- -03155 

.03889 

90 

1.00000 

1 .00000 

.15665 

•07833 

- .03442 

-.03442 

.04240 



I 
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TABLE 8.- RADIAL DEFLECTIONS v/R AT VARIOUS VAIDES OF 9 


[Family II. k' = 0 . 6 } t' = 2 .5] 


1 

2 

3 

1 + 

5 

6 

7 

0 

(deg) 

cos 20 

cos 1+9 

-0.07166 X (2) 

-0.00274 x (3) 

0.02503 

I - ( 4 ) + ( 5 ) + ( 6 ) 

0 

1.00000 

1.00000 

-0 .07166 

-0 .00274 

0 .02503 

- 0 .04937 

10 

.93969 

.7660U 

- .06734 

- .00210 

02503 

- .04441 

20 

.76601+ 

.17365 

- .051+89 

- .00048 

02503 

-.03034 

30 

.50000 

-.50000 

-.03583 

.00137 

02503 

- .01217 

1+0 

.17365 

-.93969 

-. 0121 + 1 + 

.00257 

.02503 

.01516 

50 

-.17365 

-.93969 

. 0121 + 1 + 

.00257 

.02503 

.o4oo4 

60 

-.50000 

-.50000 

.03583 

.00137 

.02503 

.06223 

70 

- .76601+ 

.17365 

. 051+89 

-.00048 

.02503 

.07944 

80 

-.93969 

. 76601 + 

. 0673 ^ 

- .00210 

.02503 

.09027 

90 

-1.00000 

1.00000 

.07166 

-.00274 

.02503 

•09395 


ON 

VO 
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TABLE 9-- ANGULAR DISPLACEMENTS 0 AT VARIOUS VALDES OP 9 


|Family II* k' = 0*6; t 



1 

2 

3 

4 

5 

6 

7 

e 

(deg) 

sin 20 

sin 40 

0.02751 x (2) 

-0.00809 X (3) 

0 = (4) + (5) 

(rad..) 

0 

(deg) 

0 

0 

0 

0 

0 

0 

0 

10 

•34202 

.64279 

.00941 

- .00520 

.00421 

.24 

20 

.64279 

•98481 

.01768 

- .00797 

.00971 

.56 

30 

.86603 

.86603 

.02382 

- .00701 

.01681 

.96 

40 

•98481 

•34202 

.02709 

- .00277 

.02432 

1.39 

50 

•98481 

-.34202 

.02709 

.00277 

.02986 

1-71 

60 

.86603 

- .86603 

.02382 

.00701 

.03083 

1.77 

70 

.64279 

-.98481 

.01768 

.00797 

.02565 

1.47 

80 

•34202 

-.64279 

.00941 

.00520 

.01461 

.84 

90 

i 

0 

0 

0 

0 

0 

0 
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TABLE 10.- BENDING -MOMENT DISTRIBUTION 


o 
1 > 


Family II- k ' = 0 .6; t ’ = 2 -5j § ^ = 1-8464 

K 


1 

2 

3 


5 

6 

7 

8 




/ F 'Y" J 1 

*0 l 

V 

v r 

M _L 

9 

r 

(zf 

W ' S e 2 

q. b 2 

E 

R R 

q. e 2 

(deg) 

E 

W 

= (3) - 1.81*61* 

= (4)/2 

from table 8 

= (2) x (6) 

= (5) + (7) - 0.01339 

0 

1.60000 

2.56000 

0 .71360 

0.35680 

-0 .01*937 

-0 .07899 

0 .26442 

10 

1.5851*5 

2.51365 

.66725 

.33363 

- .01*1*1*1 

- .07041 

.24983 

20 

1.51*252 

2-37937 

-53297 

.2661*9 

- .03031* 

- .04680 

.20630 

30 

1-1*7356 

2.17138 

.321*98 

.1621*9 

-.01217 

- -01793 

.13117 

1*0 

1.38226 

1.910 64 

.061*21* 

.03212 

.01516 

.02096 

.03969 

50 

1.27376 

1.6221*6 

-• 2239 1 * 

-.11197 

.01*001* 

.05100 

- .07436 

60 

1.151*69 

1-33331 

- -51309 

-.25655 

.06223 

.07186 

-.19808 

70 

1.06418 

1.1321*8 

- .71392 

-.35696 

.0791*1+ 

.08454 

- .28581 

80 

1.0151*2 

1.03108 

- .81532 

-.1+0766 

.09027 

.09166 

-.32939 

90 

1 .00000 

1 .00000 

- .81*61*0 

-.1*2320 

.09395 

•09395 

1 

-.34264 
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TABLE 11.- DEFLECTIONS u AND v AT 
VARIOUS VALDES OF 0 


[Family II. k’ = 0.6; t' = 2-5^ 


1 

2 

3 

0 

(deg) 

u/R 

v/R 

0 

0 

- 0 .05064 

10 

.00001 

-.04649 

20 

.00446 

-•03553 

30 

.01549 

- .02158 

40 

.03270 

- .00920 

50 

.05300 

-.00182 

60 

.07165 

0 

70 

.08424 

0 

80 

.08726 

0 

90 

.08939 

0 
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Original 

Deflected 


P-* 



► P 


Figure 1.- Initially curved bar in original and in 

deflected shape. 


y 



Figure 2.- Forces and moments 
on a pressure -loaded ring. 


y 



Figure 3. - Angular rotation and 
radial deflection of a point on 
a ring after loading. 


Figure 4.- 



Coordinate system used in beam deflection. 



y 



Figure 5.- Shapes of rings used in numerical examples. 
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a R in terms of a^/EI 
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I lil I I 1 1 I I 

0 .1 .2 .3 .4 


Figure 6 .- Fourier coefficients a n 

9 

parameter k. k = 


for various values of 



in terms of a" 5 /El 
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Figure 7 .- Fourier coefficients b R for various values of 
parameter k. = 1 - . 


L 


in terms of a 4 /EI 
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Figure 8.- 


Fourier coefficients c n 
parameter k. = 1 


for various values of 



in terms of 
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0 .1 .2 .3 .4 

Figure 9.- Fourier coefficients o for various values, of 

n 

parameter k. k z = (r) • 


in terms of R Z /E1 
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1U 



parameter k' . k' = H/R. 


in terms of Rr/EI 


80 
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k ' 

Figure 11.- Fourier coefficients b n for various values of 
parameter k T . k' = H/R. 


in terms of R 4 /EI 
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Figure 12 .- Fourier coefficients c n for various values of 
parameter k ’ . k ' = H/R. 


in terms of 1/R 
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0 .2 .4 .6 

k’ 


Figure 13.- Fourier coefficients e n for various values of 
parameter k'. k’ = H/R. 
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Figure 15.- Values of J and S. 


Family n. k ' 


H/R. 
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(T T7»T 1 

Figure 18 .- Values of A 4 /a for various ratios of — . t = — — . 

4/ El q a 3 


k 


2 


= 1 



2 
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k' = H/R. 
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k' 

Q El 1 

Figure 20.- Values of Ag/R for various ratios of t' = — 


k' = H/R. 


90 
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0 .2 .4 .6 


k* 


Figure 21.- 


Values of A^/R for various ratios of -3- . 


t' = -L . k' = H/R. 

q r 3 



L 
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0 10 20 30 40 50 60 " ' ’ 70 80 90 

9 , deg 

Figure 23.- Values of w/a for various positions around ring from t = 0.20 to t = 0.50. 
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E=a 


CO 

> 



Figure 26.- 


0, deg 

Values of w/R for various positions around ring from t * = 1.5 to t'=4.0. 


t* 


Si±. k . 

q r 3 ’ 





VO 

VJ1 
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ON 


w/R 



Figure 27.- Values of w/R for various positions around ring from t'=1.5 to t' = 4.0. 
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20 


30 


60 


70 


80 


40 50 

9, deg 

Figure 29.- Values of 0 for various positions around ring from t = 0.20 to t = 0.50. 

t = f 73 ’ k2 = 1 ‘(I) 2 - ai - 
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Figure 44.- General view of test setup. 
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Figure 45.- Schematic drawing of test setup. 
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Figure 48.- Pressure bag used to load ring. (Uniform load could 

not be obtained.) 
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by present method, a = 27.74 in.; b = 23.24 in.; El = 105,410; q = 10.12 lb/in.; 

k 2 = 1 - (|) 2 = 0.296; t = \ = 0.488. 
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by present method, a = 27.74 in.; b = 23.24 in.; El *= 105,410; q = 10.12 lb/in.; 

k 2 = 1 - ( - = 0.296; t = — i = 0.488. 
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Figure 52.- Comparison of experimental values of w/R and values of w/R computed 
by present method. H = 4.77 in.; R = 23.85 in.; El = 105,410; q = 10 lb/in.; 

k ' = 3 = 0.2; t’ = ~ — = 0.794. 
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Figure 53.- Comparison of experimental values of 0 and values of 0 computed 
by present method. H = 4.77 in.; R = 23.85 in.; El = 105,410; q = 10 lb/in.; 

k' = | = 0.2; t»= * = 0.794. 
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by present method. H = 4.77 in.; R = 23.85 in.; El = 105,410; q = 10 lb/in.; 

k ' = S * 0.2; t' = — -L = 0.794. 
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Figure 55.- 



Increase of arc length due to angular displacements. 



Figure 56.- Increase of arc length due to radial deflections. 
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